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Differential elliptic flow (112) for <fi mesons and (anti)deuterons (d)d is measured for Au+Au colli- 
sions at -yijvjv = 200 GeV. The V2 for <j> mesons follows the trend of lighter 7T and K mesons, 
suggesting that ordinary hadrons interacting with standard hadronic cross sections are not the pri- 
mary driver for elliptic flow development. The V2 values for (d)d suggest that elliptic flow is additive 
for composite particles. This further validation of the universal scaling of vi per constituent quark 
for baryons and mesons suggests that partonic collectivity dominates the transverse expansion dy- 
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An important goal of current ultra-relativistic heavy 
ion research is to map out the accessible regions of the 
Quantum Chromodynamics (QCD) phase diagram. Cen- 
tral to this goal, is the creation and study of a new 
phase of nuclear matter - the Quark Gluon Plasma 
(QGP). Thermalization and de-confinement are impor- 
tant properties of this matter, believed to be produced 
in heavy ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) P,S,S]. 

Detailed elliptic flow measurements provide indispens- 
able information about this high energy density mat- 
ter [1, H, @, 0, S|- Such measurements are character- 
ized by the magnitude of the second-harmonic coeffi- 
cient V2 — ^e l2 ' ¥>p_ * i?p )), of the Fourier expansion of 
the azimuthal distribution of emitted particles. Here, 
tf p represents the azimuthal emission angle of a particle, 
$AP is the azimuthal angle of the reaction plane and the 
brackets denote statistical averaging over particles and 
events (sllio). 

At RHIC energies, there is now significant evidence 
that elliptic flow, in non-central collisions, results from 
hydrodynamic pressure gradients developed in a locally 
thermalized "almond-shaped" collision zone. That is, 
the initial transverse coordinate-space anisotropy of this 
zone is converted, via particle interactions, into an az- 
imuthal momentum-space anisotropy. Indeed, when plot- 
ted as a function of the transverse kinetic energy KEt = 
rriT — m divided by the number of valence quarks n qi 
of a given hadron (n q = 2 for mesons and n q = 3 
for baryons), V2/n q shows universal scaling for a broad 



range of particle species [111, QJl, [13| (rriT is the trans 



verse mass). This has been interpreted as evidence that 
hydrodynamic expansion of the QGP occurs during a 
phase characterized by (i) a rather low viscosity to en- 
tropy ratio 77/s @, H, 0, [3, 15 1 and (h) independent 
quasi-particles which exhibit the quantum numbers of 
quarks 0, EH E3, E3, EI 0]. A consensus on the 
detailed dynamical evolution of the QGP has not been 
reached [a. Il5j|. 

Elliptic flow measurements for heavy, strange and 
multi-strange hadrons 21, |22| can lend unique insight 



on reaction dynamics. Here, we use differential vi mea- 
surements for the 4> meson and the deuteron to address 
the important question of how the existence of a hadronic 
phase affects V2, i.e whether or not elliptic flow develop- 
ment is dominantly pre- or post-hadronization. 

The (j> meson is comprised of a strange (s) and an 
anti-strange (s) quark and its interaction with hadrons 
is suppressed according to the Okubo-Zweig-Izuka (OZI) 
rules [23j. One consequence of this is that the <j) meson 
is expected to have a rather small hadronic cross section 
with non-strange hadrons (~ 9 mb) [3, [2^, [26[ . Such 
a cross-section leads to a relatively large mean free path 
A^, when compared to the transverse size of the emitting 
system IH 2Q|. Thus, if elliptic flow was established in a 
phase involving hadrons interacting with their standard 



hadronic cross sections (post-hadronization), one would 
expect V2 for the <f> meson to be significantly smaller than 
that for other hadrons (e.g. p and ir). If V2 is established 
in the phase prior to hadronization, the 4> meson pro- 
vides an important benchmark test for universal scaling 
in that its mass is similar to that of the proton and the 
A baryon, but its V2 should be additive with respect to 



the i>2 of its two constituent quarks (i.e. 



2). There- 



fore, a detailed comparison of the V2 values for the 4> 
meson with those for other particle species, comprised of 
the lighter u and d quarks or the heavier charm quark 
c, can provide unique insight on whether or not partonic 
collectivity plays a central role in reaction dynamics at 
RHIC [H,S,[2|. 

The deuteron is a very shallow composite p + n bound 
state, whose binding energy (~ 2.24 MeV) is much less 
than the hadronization temperature. Thus, it is likely 
that it would suffer from medium induced breakup in the 
hadronic phase, even if it was produced at hadronization. 
In fact, recent investigations [H [3(| suggest that (pn)pn 
coalescence dominates the (anti) deuteron (d)d yield in 
Au+Au collisions. Thus, V2 measurements for (d)d also 
provide an important test for the universal scaling of el- 
liptic flow [27[ in that its V2 should be additive; first, with 
respect to the V2 of its constituent hadrons and second, 
with respect to the V2 of the constituent quarks of these 
hadrons, i.e. n q — 2 x 3. 

In the 2004 running period the PHENIX detector [31 1 
recorded « 6.5 x 10 s minimum-bias events for Au+Au 
collisions at ^/snn — 200 GeV. The collision vertex 
z (along beam axis) was constrained to \z\ < 30 cm 
of the nominal crossing point. The event centrality 
was determined via cuts in the space of Beam-Beam 
Counter (BBC) charge versus Zero Degree Calorimeter 
energy [32|]. In the central rapidity region (|?7| < 0.35) 
the drift chambers, each with an azimuthal coverage 
Aip = 7r/2, and two layers of multi-wire proportional 
chambers with pad readout (PCI and PC3) were used for 
charged particle tracking and momentum reconstruction. 
The time-of-flight (TOF) and lead scintillator (PbSc) de- 
tectors were used for charged particle identification [g, 0] ■ 

Time-of-flight measurements from the TOF and PbSc 
were used in conjunction with the measured momen- 
tum and flight-path length, to generate a mass-squared 
(m 2 ) distribution [33j ]. A track confirmation hit within 
a 2.5(7 matching window in PC3 or TOF/PbSc served 
to eliminate most albedo, conversions, and resonance 
decays. A momentum dependent ±2(7 cut about each 
peak in the m 2 distribution was used to identify pions 
(7r ± ), kaons (K^), (anti)protons ((p)p), and {(d)d) in 
the range 0.2 < p T < 2.5 GeV/c, 0.3 < p T < 2.5 GeV/c, 
0.5 < p T < 4.5 GeV/c, and 1.1 < p T < 4.5 GeV/c re- 
spectively in the TOF, and to identify in the range 
0.3 <p T < 1-5 GeV/c in the PbSc. This gives ~ 59000 
d + d. An invariant mass analysis of the </> — > K + K~ 
decay channel yielded ~ 340000 <j> mesons with relatively 
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good signal to background (14 - 42% for the mass win- 
dow | mint, | = 5 MeV/c 2 about the <\> meson peak) over 
the range 1.0 < pt < 5.5 GeV/c for K + K~ pairs. 

The reaction plane method [6| was used to correlate the 
azimuthal angles of charged tracks with the azimuth of 
the event plane $2, determined via hits in the two BBCs 
covering the pseudo-rapidity range 3 < \r]\ < 3.9. The 
large rapidity gap A77 > 2.75 between the central arms 
and the particles used for reaction plane determination 
reduces the influence of possible non-flow contributions, 
especially those from di-jets [341 ]. 

Charge averaged values of v-i = 
(cos(2(^ p - $ 2 ))>/(cos(2($ 2 - ®rp))) were evalu- 
ated for tt^ 1 j K , (p)p, and (d)d. Here, the denominator 
represents a resolution factor which corrects for the 
difference between the estimated $2 and the true az- 
imuth $>rp of the reaction plane @, [35| . The estimated 
resolution factor of the combined reaction plane from 
both BBC's has an average of 0.33 over centrality, with a 
maximum of about 0.42 in mid-central collisions [il fl2T|. 
The associated systematic error is estimated to be ~ 5% 
for 7T , K , and (p)p. A px dependent correction 
factor (~ 5 — 11%) was applied to the v 2 values for 
(d)d, to account for background contributions to the 
(anti)deuteron peak (signal) in the m 2 distributions (see 
dashed-dot curve in Fig.[TJi): 



(d)d, , 
v 2 (p T ) 



= {vl +ba ( PT )-(l-R)v b 2 \p T ))/R, (1) 

where «2 +bs '(pt) is the measured t>2 for (d)d + 
background at a given px, R is the ratio sig- 
nal/ (signal+background) at that pt, and v\ g (pr) is the 
elliptic flow of the background evaluated for m 2 values 
outside of the (d)d peaks. 

Extraction of the elliptic flow values for the <j> me- 
son («2 ) followed the invariant mass (rrii nv ) method [36| . 
For each event, m,- n „, p^ lr , and tpP air for each K + K~ 
pair were evaluated. Then, for each p^? lr bin, v 2 air = 
(cos (2(tpP mr - $ 2 ))> was evaluated as a function of m,™ 
as shown in Fig.[T}3. The value v% (pt) was then obtained 
from ^2 air {rriinv ) via an expression similar to Eq. [TJ 

v% ml (m mv ) = v$ R(m inv )+v 2 s (m inv )(l- R(m inv )), (2) 

where R(m mv ) = N$(mi nv ) /[N^[m inv ) + N bg (m inv )} 
and N$(mi nv ) and Nb g (m inv ) are distributions for the 
<f> meson and the combinatoric background, respectively. 
N^irriinv) is obtained from the distribution N pa i r (mi nv ) 
of K + K~ pairs from the same event (foreground); 
Nbg(mi nv ) is the distribution of pairs obtained from dif- 
ferent events with similar centrality, vertex, and event 
plane orientation [37J. FigureQJa) shows a representative 
example of the latter distributions for 1.6 < p^ lr < 2.7 
GeV/c and reaction centrality 20-60%. A clear peak sig- 
naling the <j) meson is apparent in the foreground distri- 
bution for rriinv ~ 1.02 GeV/c 2 . The background dis- 
tribution was normalized to that for the foreground in 
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FIG. 1: (a) m inv distributions for foreground (points) 
and background (dashed-line) K + K~ pairs (py° lr = 1.6 — 
2.7 GeV/c) for 20-60% central Au+Au collisions. (b) 
rriinv distribution after subtraction of the background; (c) 
(cos {2(ip pal ' r — $2))) vs. m,i nv ; the solid line is a fit to the 
data with Eq.[2] (d) m 2 distribution for d, d for pr = 1.6 — 2.9 
GeV/c. 



the region 1.04 < rriinv < 1-2 GeV/c 2 and subtracted to 
obtain the N^(mi nv ) distribution shown in Fig. [l^b); a 
relatively narrow <p meson peak is apparent. 
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FIG. 2: (color online) Comparison of differential vafjJr) for <j> 
mesons, (d)d, tt , K , and (p)p (as indicated). Results are 
shown for 20-60% central Au+Au collisions. 

Determination of the ratio R(mi nv ) was facilitated by 
fitting this distribution with a Breit-Wigner plus a linear 
function, as shown by the solid curve in Fig.[ljc). To en- 
sure robust V2 extraction, the combinatorial background 
was constructed such that v \f mix ^{minv) gave the same 
value as v^" (rrii nv ) for m,„„ values not associated with 
the 4> meson peak. Values for v% were extracted via direct 
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FIG. 3: (color online) (a) vi vs KEt for several identified par- 
ticle species obtained in mid-central (20-60%) Au+Au colli- 
sions, (b) V2/n q vs KEr/riq for the same particle species 
shown in panel (a). The shaded bands indicate systematic 
error estimates for (d)d and <j> mesons (see text). 



fits to the v^" {mi nv ) distribution for eachp^" r selection 
(cf. Eq.[2|). That is, v^ 6 {mi nv ) was parametrized as a lin- 
ear or quadratic function of to.;„„ (depending on the p T a %r 
bin) and was taken as a fit parameter. 

The accuracy of the extraction procedure was verified 
by checking that the m inv dependence of the sine co- 
efficients, v p ™ r (m mv ) = {sin (2{ipP air - $ 2 ))), were all 
zero within statistical errors. An alternative "subtraction 
method" [H, [39| , in which the raw <p meson yield distri- 
bution dN/d(ip,p — $ 2 ) was extracted and fitted with the 
function N(l + Iv^cos (2(<^ — $2))), also showed good 
agreement, albeit with larger error bars; N is an arbitrary 
normalization constant. 

The differential V2(pt) obtained for (a)d and the me- 
son, for centrality 20 - 60%, are compared to those for 7r ± , 
if ± , and (p)p in Fig. [2j This centrality selection was so 
chosen to (i) maximize the <p meson signal to background 
ratio over the full range of pr bins and (ii) enhance the 
distinction between baryon and meson w 2 in the interme- 
diate pr range. The shaded bands for (d)d and the (j> 
meson indicate systematic errors (~ 6 — 15%), primarily 
associated with the determination of R and R(rrii nv ), v^ 9 
and V2 9 (minv) (cf. Eqs. [T]and[2]), and fitting. 

The values for shown in Fig. [2] are as much as a 
factor w 2.5 lower than those for 7r ± at Iowpt- This mass 
ordering pattern reflects the detailed expansion dynam- 
ics of the created matter. As a first test of whether or 
not V2 for (d)d is additive with respect to its constituent 

hadrons, vf^* 3, j2 vs. pt/2 is compared to v% p vs. pr- 
Within errors, they show good agreement as would be 

expected if is additive. Another salient feature of 
the results shown is the large magnitude of v-i for the <p 
meson, which gives an initial indication that significant 
flow development occurs prior to hadronization. 

The left and right panels of Fig. [3] compare the un- 



sealed and scaled results (respectively) for V2 vs. KEt 
for 7r ± , K^, (p)p, (d)d, and the (f> meson, in 20-60% cen- 
tral Au+Au collisions. The left panel clearly shows that, 
despite its mass which is comparable to that for the pro- 
ton, V2(KEt) for the 4> meson follows the flow pattern 
of the other lighter mesons (ir and A"), whose cross sec- 
tions are not OZI suppressed. A similar pattern is also 
observed for the V2(KEt) values inferred for D mesons 
(comprised of charmed quarks) from non-photonic elec- 
tron measurements We interpret these observa- 
tions as an indication that, when elliptic flow develops, 
the constituents of the flowing medium are not ordinary 
hadrons interacting with their standard hadronic cross 
sections. Instead, they may indicate a state in which 
partonic collectivity dominates the transverse expansion 
dynamics of light, strange, and charmed quarks via a 
common velocity field. 

Interestingly, the V2(KEt) results shown for the (d)d 
and the 4> meson are essentially identical at low KEt 
{KEt ^ f GeV), and are in good agreement with those 
for other charged hadrons, including the pion with a mass 
^ f 3 times smaller than the deuteron. This strengthens 
the earlier finding that, for low KEt, all particle species 
exhibit the same V2 irrespective of their mass [Tf|, EH, 
UU. The expected difference between (d)d and (p)p for 
KEt ^ f GeV is not tested in Fig. [31 due to the limited 
KEt range of the (d)d data. 

The right panel of Fig. [3jshows the results for a valida- 
tion test of universal scaling for V2(KEt) of baryons and 



mesons 



IJ,ll2j,lf3|. The value n q = 2 x 3 is used for [d)d 
to account for its composite (p + n) nature. The scaled 
results shown in Fig. [3Jd clearly serve as further vali- 
dation for the experimentally observed universal scaling 
of V2 for baryons and mesons [O, HH, [T3|. This find- 
ing lends strong support to the notion that the high en- 
ergy density matter, created in RHIC collisions, comprise 
a pre-hadronization state that contains the prerequisite 
quantum numbers of the hadrons to be formed. Thus, 
it appears that partonic collectivity dominates the ex- 
pansion dynamics of these collisions. The special role of 
KEt as a scaling variable is under investigation. 

In summary, we have presented differential V2 measure- 
ments for the cj) meson and deuteron, and have compared 
them to those for other mesons and baryons. For a broad 
range of KEt values, the differential V2{KEt) for the 4> 
meson follows the flow pattern for other light mesons 
whose cross sections are not OZI suppressed. The com- 
posites (d)d follow the flow pattern for baryons with V2 
values which are additive. When V2/n q is plotted as a 
function of the transverse kinetic energy scaled by the 
number of valence quarks (ie. KEt /n q ), universal scal- 
ing results for all particle species measured. These obser- 
vations suggest that the transverse expansion dynamics 
leading to elliptic flow development cannot be understood 
in terms of ordinary hadrons interacting with their stan- 
dard hadronic cross sections, but rather in terms of a 
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pre-hadronization state in which the flowing medium re- 
flects quark degrees of freedom. 
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